Oscillating Superfluidity of Bosons in Optical Lattices 
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Following a suggestion by Orzel et al [1], we analyze bosons in an optical lattice undergoing a sudden 
parameter change from the Mott to superfluid phase. We introduce a modified coherent states 
path integral to describe both phases. The saddle point theory yields collective oscillations of the 
uniform superfluid order parameter. We calculate its damping rate by phason pair emission. In two 
dimensions the overdamped region largely overlaps with the quantum critical region. Measurements 
of critical dynamics on the Mott side are proposed. 
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With recent experimental developments of ultra cold 
atoms in optical lattices, the fascinating phenomena of 
Bose-Einstein Condensation have entered the domain of 
strong interactions [1,2]. Macroscopic quantum states 
can be effectively manipulated and time evolution of or- 
der parameters (OP), adiabatic [2] or non-adiabatic [1], 
can be probed by varying the optical lattice parameters. 

In one such experiment, the strength of a three di- 
mensional optical lattice potential was tuned to induce a 
quantum phase transition between a Mott insulator and 
a superfluid of bosons [2] . This phase transition has been 
extensively analyzed theoretically [4-8] and numerically 
[9] . The two phases are characterized by markedly differ- 
ent many body states. The Mott phase, at large lattice 
potential barriers, is well described by definite real space 
occupation numbers. The compressible superfluid phase, 
on the other hand, sustains long range phase order. This 
phase is detected by self interference patterns after the 
gas is released from the trap. 

In an interesting proposal, Orzel et al [1] suggested 
the possibility of observing OP time evolution. Basi- 
cally, the bosons are prepared in the number squeezed 
Mott state, and then the potential is suddenly reduced 
into the superfluid phase. The consequent evolution of 
the superfluid order can be deduced from the intensity of 
interference patterns appearing when the atoms are re- 
leased from the trap at sequential times. This would open 
up exploration of a new regime of macroscopic quantum 
dynamics [10-12]. The initial questions which come to 
mind are: (i) Could coherence oscillations be observed? 
(ii) What would be the time scale of superfluid OP evolu- 
tion? (iii) What would be the damping mechanism, and 
damping rate of such effects? 

These are the primary issues addressed in this paper. 
We derive the effective Hamiltonian of the superfluid OP 
starting from interacting bosons in a periodic lattice. We 
find a variational bosonic representation which describes 
the phase diagram and treats the elementary excitations 
on both the Mott and superfluid phases. In the Mott 
phase, the two degenerate gapped excitations become 



gapless at the transition. The superfluid phase reduces 
to a relativistic Gross-Pitaevskii action, with one (gap- 
less) phase mode and one (gapped) amplitude mode. We 
obtain a semiclassical solution of a macroscopically oscil- 
lating superfluid order and calculate its relative damp- 
ing rate. This provides an estimate of the experimental 
regime where such oscillations should be visible. In two 
dimensions this region largely overlaps with the quantum 
critical region, as we estimate from Ginzburg's criterion. 
We end by commenting on critical dynamics, and how 
they might be observed. 

The Bose Hubbard model (BHM) describes interacting 
bosons in an optical lattice, 
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where n, is the boson occupation on site i. The tun- 
neling J and interaction U are known functions of the 
microscopic forces [3]. At integer fillings, n = 1,2,..., 
the BHM exhibits quantum phase transitions. For large 
tunneling (weak optical potential barriers) Jfi > > U the 
ground state is a superfluid (Bose Condensed) with long 
range phase order. Below a critical tunneling strength 
J < J c (n), bosons are localized in incompressible (inte- 
ger occupations) Mott phases. 

In the vicinity of the Mott phase, number fluctuations 
are small. An effective Hamiltonian truncated into the 
subspace of lowest local number states \n— 1), |n), |n+l), 
captures the essential correlations around the transition. 
The reduced Hilbert space can be represented by three 
commuting t-bosons |fl+a) = i^jO), a — 1,0, — 1, which 
obey the holonomic constraint ^ Q t ai t a i = 1. In this 
subspace, the bosons of the BHM are represented by a\ — 

\fh~ 4i*-i« + V™ + 1 4i*Q»- 

At large n [13], the effective Hamiltonian assumes a 
particularly simple pseudospin-one form 
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St = V2(4*-h + tltoi) S? = t\ ti t lti - tl hi t- lti . (2) 

The complex supcrfluid OP field breaks planar spin 
symmetry ^(xj) = ^/E(St)- It is tempting to describe 
the action of i? e // using spin one coherent states [14]. 
However, the Mott ground state is perturbatively con- 
nected to the 0(2) rotationally invariant state JJ { \Sf — 
0) . Thus it is difficult to describe this phase as a saddle 
point of a spin coherent states path integral. 

Alternatively we can use modified coherent states de- 
fined by 

\n(e, V ,cj>,x)) - Us(0/2)4 + e"'sin(0/2) 



EL %% I ^ ) • The remaining b* a and 6^ bosons create fluc- 
tuations about this state. The next step is to apply the 
constraint in (6), and eliminate bo from the Hamiltonian 
[15]: 
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x (e^ sin(x/2)4 + e - ** cos(x/2)iO 1 | ) . (3) operators b a and 6, 



The mean field theory is similar to previous variational 
approaches [5,6]. A homogeneous variational wave func- 
tion which captures both phases is | $ m / ) = J}, I ) ■ 
= describes the Mott phase, while 9 > has a super- 
fluid OP oc ^/EsinO. The variational energy per site 
is 
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x(l + n 1 sin 2 (x/2) + \J \ + n^ 1 sin* cos 277) (4) 

where z is the lattice coordination number. The Mott 
phase boundaries found by minimizing e var are given by 



^=-8^ ± W i -U i+ i 



(4nu)~ 2 , (5) 



where u = U / (AJnz). For large n (and commensurate 
density) the Mott transition occurs at fi — and u = 1. 
For n = 1 it occurs at £/ = 5.8z J. 

The error incurred by truncation to three states per 
site is estimated by comparing (3) with a variational an- 
zats which indues 11 occupation states. Even for n » 1 
we find a wide regime (u > 0.4) in which the probability 
weight of states out side the truncated Hilbert space is 
less than 1%. The particle number fluctuation Sn < 0.6 
in this range and the error in Sn is less than 10%. 

To keep the presentation simple, we focus on the limit 
of large occupation numbers. This amounts to treating 
(2), which we believe retains the correct qualitative dy- 
namics even for low occupations. Excitations. Consider 
the canonical transformation 

b 0l = cos(6»/2)t + sin(0/2)(t u + t- U )/V2, 
b m = sm(<9/2)i - cos(0/2)(t u + t- U )/s/2, 
b vi = (t u - t_ii)/A (6) 

with local constraints ^2 m b mi b m i = 1. The mean field 
variational state is simply the Fock state | <J> m / ) = 



Truncation at quadratic order is valid provided 
(btxba), (btb v ) <C 1, which can be tested self consistently 
in the regime of interest. 

Now we can expand Eq. (2) in terms of the fluctuation 
to obtain a harmonic hamiltonian 
with normal and anomalous terms. It is diagonalized by 
a standard Bogoliubov transformation to obtain 



Hfiuc = ^2 w ™k/3^ k /3 m ,k to = 1, 2. 
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In the superfluid phase, there is an amplitude mode and 
phasons 



w a (k) = 2z Jny/l - u 2 7 k « V c 2 k 2 + A 2 , 
w v (k) = Jzn(l + u)y/l - 7k, « cjkj, 



(9) 



where 7 k = l/zE< e 
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2 j2;n-\/l — u 2 . The massive amplitude mode u a softens 
at the Mott transition and becomes degenerate with u) v . 
In the Mott phase, the two modes 
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are gapped, representing particle and hole excitations. 
The local density of fluctuations are found to be rela- 
tively small (b mi b m i) < 0.08 even at the critical point. 
This measures the accuracy of the variational state (3) at 
least for the short wavelength correlations. In the critical 
region, which will be estimated later, we do not expect 
this mean field theory to yield correct exponents for the 
vanishing of the gap or the long wavelength correlations. 

Using the modified coherent states | O ) , we can con- 
struct a path integral for the evolution operator as fol- 
lows. A resolution of identity is found to be of the form: 
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with the invariant measure A4(9) = C cos 9(3 cos 8 — 1). 
It is also straightforward to calculate the kinetic term 

=isin 2 (6>/2)(?7-cosx^) = -«r(t). (12) 

(11) and (12) are the necessary ingredients for the path 
integral 
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U(t) = J TXlM{9)ex${% J dt'[r{t')-H{t')]} (13) 

We now focus on the commensurate case fi = 0. An 
action in terms of ^ = ^/Esm9e~ ^^p is obtained by in- 
tegrating over the massive fields n and x- Expanding to 
fourth order in \*5f\ and taking the continuum limit we 
arrive at a "relativistic" Gross-Pitacvskii action 
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This derivation is valid to second order in the dimension- 
less distance to the critical point |1 — u\. 

Saddle point collective oscillations. The time evolu- 
tion of the system after a sudden change of parame- 
ters from the Mott to the supcrfiuid state is determined 
by the saddle point of the action (14). The equation 
of motion for the OP, rescaled by its equilibrium value, 
I *o I = ^/2n(l — u)/u is given by 



* = c 2 V 2 * + -A 2 *(l-|*| 2 ). 



(15) 



The constants A = 2\/2 J 'zfiy/ T — u and c = 2Jn v / i are 
identical to the expressions below (9) to leading order in 
1 — u\. For a uniform field configuration (15) is readily 
integrated to give the non linear oscillations shown in 
Fig. lb. whose time scale is set by A. Fig. la shows the 
motion of the OP on the potential landscape. 
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FIG. 1. Oscillations of the superfluid order parameter vp 
in a system prepared in the Mott state, (a) Mean field en- 
ergy, where the equilibrium ground state is marked by point 
(2). (b) Solutions of the saddle point Eq. (15), which ex- 
hibits macroscopic oscillations. The dashed line indicates the 
superfluid OP in the ground state (2) . 



Restriction to a uniform field is not justified apriori. 
In particular, topological defects may be trapped by the 
so called Kibble-Zurek mechanism [16]. However, the 
number of trapped vortices can be diminished. If the 
initial Mott state is close to the transition, a large cor- 
relation length £m determines the distance between seed 
vortices. If on the other hand we start deep in the Mott 
phase, we should consider initially zero OP with small 
random fluctuations, uncorrelated on the scale of a lat- 
tice constant. This amounts to a seed vortex on almost 
every plaquette. However, we argue that only few will 
survive the first stage of evolution. To describe the ini- 
tial growth of the OP consider the linearized version of 
(15) around 'J = 0, whose eigenmodes are easily found to 

be w(k) = J (cfc) 2 — |A 2 . Since it is the fastest growing, 



the uniform k = mode will dominate the development 
of an OP. More over, fluctuations with k > l/(£\/2) do 
not grow at all. This implies that defects with the same 
topological charge must be separated by at least £ = c/ A, 
a large distance in the regime of interest, not far from the 
transition point. The impact of the few remaining vor- 
tices on the evolution is an interesting open issue which 
deserves further study. 

Damping of the oscillations. The collective oscillations 
in Fig. 1 are in fact a macroscopic occupation (in a co- 
herent state) of the zero wave vector amplitude mode, 
w a (k = 0). Since this mode is coupled anharmonically 
to the low energy phasons, we expect a finite damping of 
the oscillations due to phason pair emission. 

Expanding the action in Eq. (14), up to the harmonic 
and cubic interactions, we obtain 



S = _L_ J dt'd d x[a 2 - c 2 (Va) 2 - A 2 a 2 
+ (|vI/ | 2 + 2|vI/ H [<p 2 -c 2 (V<p) 2 ] }, 
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where a = \^\ — |^o| is the linearized amplitude mode 
and c = 2y / zJn. In order to compute the damping rate 
of the oscillating field a, it is convenient to recast the 
continuum theory in operator form, using the amplitude 
and phason operators of Eq. 8 



a q =o 
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Phason pair creation is dominated by the vertex 
Hint = -j= £ V k (/?i,o$, k /j£_ k + H.c.) . 



(17) 



(18) 



By Fermi's golden rule, the damping rate is 
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At wave vector |k| = A/(2c), the vertex coupling con- 
stant is given by 



(2Jnzf 
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u(l -u) 2 . 



(20) 



Consequently the relative damping rate diverges in one 
and two dimensions as 



Qd 
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Oscillations could be observable for T/A < 1 which im- 
plies u < 0.59 for d = 1 and u < 0.73 for d = 2. In three 
dimensions, the Qd is finite at the transition. 

The damping ratio (21) can also be derived from the 
one loop correction to the longitudinal susceptibility [17]. 
Higher order terms can be resummed using a large N 
or renormalization group approach to obtain expressions 
valid in the critical region [17-19]. 

We comment that emitted phasons will eventually 
thermalize, bringing the system to a new, finite temper- 
ature, equilibrium. However, we are concerned with the 
shorter time transient of the system, regardless of its ul- 
timate equilibrium fate. 

Critical phenomena. Time dependent experiments 
such as those proposed in Ref. [1] could potentially mea- 
sure quantum critical fluctuations directly. By our mean 
field theory, the collective oscillation frequency, u>i(k = 
0), vanishes at the transition as (1 — u)* according to 
Eq. (9). We can estimate the quantum critical region 
below the critical dimension d = 3, using Ginzburg's cri- 
terion for the d+1 dimensional action (14). For d = 2, 
we have obtained |1 — u\ < 0.15, which overlaps with the 
overdamped region of the critical amplitude oscillations. 

In the critical region, the mean field gap exponents as 
read from (9), and (10) should be modified. To leading 
order in e = 3 — d the gap exponent is given in Ref [19] 
by: v= i + 0.1e. 

Experimental parameters for the Boson Hubbard 
model extracted from Ref. [1], can be translated into 
n w 50, 2nh/(2Jzn) w 0.7ms. The oscillation period 
is therefore larger than this timescale by a factor of 
1/a/I - ~u i . 

In summary, we have described the dynamics of bosons 
in an optical lattice using a modified coherent states path 
integral. This affords a unified description of both the 
supcrfluid and Mott phases. A system prepared in the 
unstable Mott state, is expected to exhibit macroscopic 
oscillations of the supcrfluid OP, with damping which in- 
creases towards the transition in one and two dimensions. 
It would be very interesting to investigate supcrfluid os- 
cillations on the Mott side of the transition, where there 
are no low energy phason modes to cause damping. Close 
to the transition, this would provide direct measurement 
of the dynamical critical exponents. 



After completing this paper we received a preprint by 
Polkovnikov etal addressing a similar Mott to SF transi- 
tion, using a dynamical rotator representation [20] . 
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